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Abstract 

We study the 1-dimensional Heisenberg antiferromagnet with s — ^ us- 
ing a Majorana representation of the s — | spins. A simple Hartree-Fock 
approximation of the resulting model gives a bilinear fermionic description 
of the model. This description is rotationally invariant and gives power- 
law correlations in the "ground state" in a natural fashion. The excitations 
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are a two-parameter family of particles, which are spin-1 objects. These 
are contrasted to the "spinon" spectrum, and the technical aspects of the 
representation are discussed, including the problem of redundant states. 



PACS numbers: 75.10.Jm, 75.50.Ee 
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I. Introduction 



The study of various representations of spins in terms of bosonic or 
fermionic operators is an old and well studied problem, reviewed nicely, for 
example, in Ref. 0. The need for exploring various representations has 
received a further impetus from the recent interest in the Heisenberg antifer- 
romagnet, as a standard model in the Resonating Valence Bond theories ||], 
i.e., models where states with no long ranged Neel order play an important 
role. The Schwinger boson representation is of very general validity, i.e., for 
any s, but the Schwinger fermionic representation is only valid for s = \ and 
gives sf — | J2a,a> c h Tac i,v'> witri tne constraint E CT 4- C ^ = 1 B H- The con- 
straint is not very easy to deal with, except in an averaged sense. Hence one 
may look for unconstrained representations. For s = | such unconstrained 
representations can be found. The so called "drone fermion" representation 
[Pi S is one °f the possibilities, where we write sf = a\(pi, s~ = 4>idi and 
sf = a\di — |, where the a's are canonical anticommuting variables, and (pi 
is a real fermion with $ = cf) and 2 = 1. Thus is a 'drone' whose only 
'job' is to make spins at different sites commute, rather than anticommute. 
In single site problems like the Kondo problem, these are useful ||. How- 
ever, this representation violates rotation invariance, since our choice of the 
z axis was arbitrary. A fully rotation invariant scheme does exist, and can, 
for example, be derived from the above, by simply rewriting the complex 
fermion a in terms of its two real components as a oc <ft x + i<p y . This leads to 
a representation with three Majorana fields, and is studied in this paper, in 
the context of the 1-dimensional Heisenberg model. 
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The plan of this paper is as follows. In Sec. II, we discuss the Majorana 
representation and the need for enlarging the Hilbert space of states in order 
to obtain a representation of the Majorana algebra. We introduce the spin-^ 
antiferromagnetic chain and its low-lying excitations in Sec. III. In Sec. IV, 
we use the Majorana representation to study the chain within a rotationally 
invariant Hartree-Fock (H-F) approximation. Since the H-F approximation 
is not unique in the general case, we require that the susceptibility calculated 
by two methods, namely, from the energy change and from the fluctuation 
spectrum, should agree. This requirement, interestingly, rules out several 
possibilities, and leads to a particular scheme which is implemented. We 
obtain a spectrum of low-lying excitations which bears a strong resemblance 
to the one discussed in Sec. III. 

We also discuss the spin of the Majorana fermion. In Sec. V, we compute 
the dynamic structure function and susceptibility, at both zero and finite 
temperatures, and contrast these with previously known results. In Sec. VI, 
we study the response of the model to uniform and staggered magnetic fields. 
We end with some concluding remarks in Sec. VII. 

II. Majorana Representation 

At each site n, we can write the spin operators S n = a n /2 in terms of 
three Majorana operators (fi n as |7], ^ 

< = ~ i<P V n<Pn, 

< = ~ i K € , 

and < = - % <f> x n 4> y n . (1) 
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(We set Planck's constant equal to 1). The operators 0^ (with a = x,y,z) 
are hermitian and satisfy the anticommutation relations 

{ C , <p b n } = 2 5 mn 5 ab . (2) 

It is interesting to note that the relation S% = 3/4 automatically follows 
from Eqs. ([l]-|2]); one does not have to impose any additional constraints at 
each site unlike the Schwinger representation ||. There is a local Z 2 gauge 
invariance since changing the sign of <p n does not affect S n . (The Schwinger 
representation has a local U(l) gauge invariance). 

For N sites with a spin-- object at each site, the Hilbert space clearly has 
a dimension 2^. We now ask, what is the minimum possible dimension which 
will allow a representation of the form given in Eqs. (|l]-|D? The answer is 
2 N+ \ N / 2 \ where [N/2] denotes the largest integer less than or equal to N/2. 
This follows from the observation that a representation for ([l|-0) is given by 

& = < 1>n , 

where [ o a m , ^ n ] = , 

and { ip m , ip n } = 2 S mn . (3) 

The minimum dimension required for a matrix representation of the spinless 
anticommuting operators ip n is 2^/^ 0. Thus the Majorana representation 
of spin-i objects requires us to enlarge the space of states; the complete 
Hilbert space of states is given by a direct product of a 'physical' space and 
an 'unphysical' one. Now suppose that the Hamiltonian is purely a function 
of the physical operators S n ; it therefore only acts on the physical states. 
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Then the unphysical part of the Hilbert space simply factorizes out; hence 
each value of the energy will have a degeneracy of 2^ N / 2 \ 

As an explicit example, consider the case N = 2. The Majorana Hilbert 
space is 8-dimensional, where the extra factor of 2 arises from the unphysical 
space. We can denote the 8 states as fff, etc. The physical operators 
Si and S2 only act on the first and second symbols respectively. The third 
symbol, which may be f or |, denotes the unphysical space. A Hamiltonian 
of the form Si • S2 only acts on the first two symbols; hence the energy 
levels will be precisely the ones of a two-site antiferromagnet, but with an 
additional degeneracy of 2 due to the third symbol. On the other hand, the 
Majorana operators can be written in the direct product form 

0i = a <g> 1 <g> a x , 
and (p 2 = 1 <8> a <g> a v . (4) 

Hence they act on the third symbol and can therefore mix up physical and 
unphysical states. 

One might worry that thermodynamic quantities like the entropy will get 
a spurious contribution proportional to iV due to the unphysical degeneracy 
of 2^/^. On the other hand, when we make approximations like the H-F 
decomposition discussed later, the physical and unphysical states get mixed 
up in an essential way. This completely changes the energy degeneracy; in 
particular, the H-F ground state is actually unique as we will see. 

We can think of 0^ as the fundamental field in our theory. Both a% 
and ipn can be written in terms of <f%, as can be seen from Eq. (|l|) and 
ipn = —i ( t ) n ( t )V n ( t ) n respectively 
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III. Antiferromagnetic Spin-| Chain 



We will now begin our analysis of a Heisenberg antiferromagnetic chain. 
The Hamiltonian is 

H = J ]T S n ■ S n+1 , (5) 

n 

where the exchange constant J > 0. We use periodic boundary conditions 
Sjv+i = (We set the lattice spacing a = 1). The spectrum of (H) is exactly 
solvable by the Bethe ansatz; in particular, the ground state energy is given 
by E Q = (— In 2 + 1/4) N J = — 0.4431iVJ. The lowest excitations are known 
to be four-fold degenerate consisting of a triplet (S = 1) and a singlet (S = 0) 
[|lT]j . The excitation spectrum is described by a two-parameter continuum in 
the (q, ou) space, where — tt < q < ir. The lower boundary of the continuum 
is described by the des Cloiseaux-Pearson relation ]10[ 

= -y- | sin q | , (6) 

whereas the upper boundary is given by 

u u (q) = 7rJ I sin - | . (7) 

We can understand this continuum by thinking of these excitations as being 
made up of two spin-| objects ("spinons") with the dispersion ||TT| 

= sin q , (8) 

where < q < tt. A triplet (or a singlet) excitation with momentum q is 
made up of two spinons with momenta q± and q2, such that < q\ < q2 < tt, 
q = qi + q2 if < q < tt, and q = q\ + q2 — 2tt if — tt < q < 0; further, 
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uj(q) = u}(qi) +uo{q2). The two-parameter continuum arises because q% can 
vary from to q/2 if < q < 7r, and from tt + q to 7r + g/2 if — 7r < g < 0. 



IV. Hartree-Fock Treatment, Ground State And Excitations 

We will now study this system using the Majorana representation. We 
write (|5|) in terms of Majorana operators to get a quartic expression, and 
then perform a Hartree-Fock (H-F) decomposition. Thus we write: 

H = - jE( <l> X n<t>X+l<t> y n+i + cycl. perm. (x,y,z) ) 

n 

ttiX^ti) + cycl. perm. (x,y,z) ] . (9) 



In principle, the H-F can be done in three different ways; however rotational 
invariance implies that only one kind of bilinear can have a non-zero expec- 
tation value in the ground state. Namely, 

g = -i{r n 0n + i ) , (10) 

where g has the same value for a — x,y, z; we also assume it to be translation 
invariant. The value of g will be determined self-consistently. We now have 
to diagonalize the quadratic Hamiltonian 

ff = y E C Cm + 3 jNJg 2 . (11) 

* a,n 

Since 0^ is hermitian, its Fourier expansion can be defined as 

e iJ " + &„, e-«" ] , (12) 



v Jy 0<g<7r 



where 

{ b ag , b\ q , } = 5 ab 5 qq , . (13) 

A similar half zone definition of the fourier transforms is possible in higher 
dimensions as well; for example, on the square lattice, we could restrict the 
sum to q x > 0. We will work with antiperiodic boundary conditions for 0° 
and even values of N in order to eliminate modes with q equal to and 
7i. This simplifies the calculation because the momenta q and — q are then 
distinct points in the Brillouin zone extending from — tt to tt. In Eq. (|i~2l) , 
q = 2n(p- 1/2) /iV, with p = 1,2, ...,N/2. In the limit N -> oo, we get 
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H = E E w(g) b\ q b aq + 3NJ ( ^ - I ) , (14) 



where the Majorana fermions have the dispersion 

u(q) = c sin q , (15) 

with c = 2g J. The H-F ground state | ) is therefore the state annihilated 
by all the b aq . Note that it is unique unlike the exact ground state, which 
has a degeneracy of 2 N / 2 within the Majorana formalism. It is curious that 
the H-F approximation gives a unique ground state which agrees with the 
degeneracy we would have obtained without the Majorana formalism. 
We now calculate ( |TUD in the H-F ground state and obtain 

9=1- ( 16 ) 

The H-F ground state energy is therefore 

3 

E HF = NJ = - 0.3040 NJ . (17) 



n 2 
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This is greater than the exact value mentioned above; indeed, one can show 
that any H-F decomposition must give an estimate for the ground state en- 
ergy which is bounded below by the exact value E a . The argument goes as 
follows. In Sec. II, we have shown that the exact ground state energy within 
the Majorana formalism is equal to the exact ground state E Q without the 
Majorana formalism, since the Hamiltonian H only acts on physical states. 
Let us therefore prove the upper bound result in the Majorana Hilbert space 
which includes both physical and unphysical states. Now the H-F calcula- 
tion is equivalent to self-consistently finding an ansatz ground state |0) and 
caclulating the expectation value of H in that. (One can show that |0) is 
an eigenstate of the Majorana fermion number operator. Hence an expecta- 
tion value of the form (ABCD) is indeed given by the H-F decomposition 
(AB)(CD) -(AC)(BD) +(AD)(BC), if the operators A, B, C and D are 
all fermionic). By the variational argument, the expectation value of H in 
any state is bounded below by E Q . 

The "spinon" spectrum has the same form as in ([8]) but has a different 
coefficient c exact = tiJ/2, whereas we find c = 4J/ir from Eq. (|I6|) . Note that 



the self consistent equation Eq. (10) also leads to Eq. (16), since we have 

2N 



E «+i = 4 E sing&^k, . (18) 



n 71 q>0 



The ground state is a singlet since it is annihilated by the total spin 
Stot = T, n Sn, for instance, by 

SL = -i E { b UK ~ h l q K) ■ (19) 

0<g<7r 

We now ask: What is the spin of a Majorana fermion? From the commutation 
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relations between S and b\ ql we find that the one-fermion state b\ q | ) 
has S = 1. More specifically, the states + itf yq ) | ), b\ q | ), and 
(^Ig — I ) have S z = 1, and —1 respectively. 

A two-fermion state can therefore have 5* = 0, 1 or 2 in general. However 
the state created by S q = J2 n S^e~ iqn , where < q < n, has the form 

s* \ o ) = - i E ( - b U b l,-k ) I o ) . ( 2 o) 

0<fc<<?/2 

and can be shown to have 5 = 1. We have thus derived the two-parameter 
continuum of triplet excitations in Eqs. ((§-0), with a prefactor 4/7T instead 
of vr/2. 

Finally, we can compute the equal-time two-spin correlation function 
G n = ( I S ■ S n I ) = for n = , 

= - ^[ 1 "(- 1 ) n ] for ^M21) 

This does not agree with the correct asymptotic behavior of G n which is 
known to oscillate as (— l) n /n. In particular, the H-F static structure func- 
tion S(q) = Y,n G n e~ tqn does not diverge as q — > tt in contrast to the correct 
S(q) which has a logarithmic divergence at tt. Note that J2 n G n = 0, as 
expected for a singlet ground state. It is interesting to observe that the 
Schwinger fermion representation yields a correlation function which only 
differs from (|21| ) by a numerical factor (see the first reference in ||). 

This Hartree-Fock state is readily generalized to finite temperatures, since 
we simply need to put in thermal population factors for the occupations of 
the fermions 

(bl q b aq ) = — * . , • (22) 

q 1 + exp(pcsmg) 

11 



Hence the self consistency condition Eq. ( |T0D together with Eqs. ( |T8| ) and 
(|22"D gives us 

----- V Sinq (23) 

9 it N ^ 1 + exposing) ' 1 ' 

It is easy to see that as T — > oo we have g — *> 0, and as T — > we have 

g — > | (1 — - ^ 2 T ), i.e., a power-law correction to the zero temperature 

'bandwidth' g. 

The H-F ground state discussed above is, unfortunately, not the one with 
the lowest energy. If we allow a dimerized expectation value g n in Eq. (|10|), 
where g n can alternate in strength from bond to bond, we find that the lowest 
energy is attained for the fully dimerized state in which g n — 1 for n even 
and for n odd (or vice versa). This corresponds to a dimerized ground state 
with an energy 

E odm = -~NJ, (24) 

o 

which is substantially lower than the earlier H-F value. There is a gap equal 
to J above the dimerized ground state. (This ground state is, of course, exact 



for the case N = 2 [T2|]). The reader may wonder why we are ignoring the 
dimerized H-F state in the rest of this paper, even though it has the lowest 
H-F energy. The reason is that we know by other methods, both analytical 
and numerical, that the correct ground state of the spin-| chain is translation 
invariant and that there is no gap above it. The H-F method is, after all, only 
an approximation, and different approximations can certainly give different 
results. We should therefore pick the H-F which agrees qualitatively with 
other methods; the ground state energy is not necessarily the best criterion 
for choosing one H-F over another. Having chosen a particular H-F on the 
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basis of certain features, we of course have to check whether it reproduces 
other features equally well. We will see in Sees. V and VI that the transla- 
tion invariant H-F yields reasonable results for the structure functions and 
susceptibilites also. 

V. Dynamic Structure Function and Susceptibility 

We recall the definition of the dynamical susceptibility 

X zz (Q,t) = i9(t) ([S*_ Q (t),S* Q ]) (25) 

/+oo 
dt X zz (Q,t) exp(toot) (26) 
-oo 

= exp(-^)-exp(-^) Q \ V><V \S' Q \„> . 

(27) 

The Zeeman coupling of a spin to a magnetic field is given by gifj,BS z B, where 
gi and /i# denote the Lande g-factor and the Bohr magneton respectively. 
The physical response function (i.e. gi\iB < S z >) is x — 9i ''HbX* 2 '(Q \ u ) ■ 
In the static limit u = 0, we have the usual thermodynamic argument for 
determining the susceptibility. If we perturb the system via the coupling 
H = H — giUsB J2n C0S (Q n )Sn, then the change in the free energy is SF = 
-gfn 2 B B 2 x zz (Q,0)9 Q , where 6 Q = 1/4 if Q ^ 0,tt, and 9 = 1/2 = W . Q 
Also recall that the static correlation function is given by 

< S - qSq > ~ J-oo - 1 - exp(-/3u;) ' (28) 



3 This factor of 9 arises because for a finite Q we drop two of the four terms in second 
order perturbation theory using momentum conservation; this neglect is disallowed exactly 
at Q — 0. 7T. 
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We will now compute the response functions in the H-F approximation. We 
begin by expressing, for < Q < it, the operator Sq in terms of the Majorana 
fields in the Heisenberg picture: 

0<q<Q 

-i 2-7T — Q — g) b xq b y ^-Q-q exp -i(u q + ^-Q- q )t 



■K—Q<q<TT 

E - 

Q<q<n 



'I E T(?> Q-Q) HqKq-Q ~ b lqKq-Q\ ^V^q ~ Uq-Cl)t ■ 



(29) 



In this equation we have introduced two real phenomenological functions 
a(a,b) = a(b,a) = a(ir — a, it — b) and 7(0,6) which are, strictly speaking, 
equal to unity from the Majorana definition of the spins. These are intro- 
duced in order to facilitate the comparison of our structure function with a 
phenomenological function proposed in Ref. | I3fl . The essential point is that 



we have assumed that the time evolution is given by the bilinear in fermions, 
our Eq. fll4|). The representation for S z _q is obtained by taking hermitean 
conjugates. Note that Sq or SLq acting on the ground state generates two 
spinons. We insert it in Eq. (|2q) , carry out the contraction of the fermions 
by Wick's theorem, and use Eq. (p2|) in the form n q =< b qa b qa > and 



n q = 1 — n q to find 



X zz (Q,uj) = £ a 2 (q,Q-q) 



n g riQ_ q - n q n Q - q 



o< q <Q ' " u i + "Q-q - w - i0+ 

0^<Q U 1 + UQ-q +UJ + Z0 + 

+ 2 E 7 2 M~Q) n ^~^-^ , (30) 

Q<q<n ^ - Uq-Q -UJ-10 + 
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This is seen to be an even function of u by using q — > n + Q — g in the last 
term. Using Eq. (^8|), we deduce that 

G ZZ {Q) =< S Z _ Q S Z Q > = Yl « 2 (^' Q ~ Q)[n q nQ- q + n q n Q _ q ] 

0<q<Q 

+ 2 7 2 (<?,9 -Q)n q - Q n q . (31) 

Q<q<n 

Let us note that at zero temperature, if we set a = 7 = 1, we get G ZZ {Q) = 
N\Q\/27r and hence the correlation function quoted in Eq. (|2T|). At the other 
extreme limit T — > 00, we replace n — n = 1/2 and find G ZZ (Q) = N/A. At 
am/ temperature, the relation n q + n q = 1 allows us to show that the sum 
rule < S*S* >= 1/4 is satisfied. 

At zero temperature, we have the static susceptibility 

X ZZ (Q,0) = 2 £ a \ Q - q) (32) 

which, in the standard situation a — 1, can be evaluated in the closed form 

y»fO 0) = - log f C ° S (?r - g)/4 1 (33) 

X W ' Uj 7rcsin(Q/2) 10g [ cos (vr + Q)/4 J " ^ 

The uniform value is 

X zz (0,0) = - = l - (34) 

TTC 4j 

The neutron scattering function which is of particular interest is found at 
zero temperature as 

SmX zz (Q,u) = 7T " 2 (<?,Q-<?) 8 {u q + u)Q- q -u). (35) 

0<g<<2 

for u) > 0. We can evaluate it in terms of the dimensionless energies u = uj/c, 

u> = 2sin(Q/2) and u < = sinQ, as 

N a 2 (a* O - a*) 

3m x»(Q, U ) = — ig-^g q ) 6( U> - u) 6(u - U< ) (36) 

c I cos(g*) — cos(<y — q*)\ 
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where g* is the solution of sing* + sin(Q — g*) = u which equals Q/2 at 



u = u > . With this we find 



sing* = - [ u — cot(Q/2)^Jv\ — u 2 
1 



cosg* = - [ wcot(Q/2) + - u 2 } . (37) 



This implies that | cos(g*) — cos(Q — g*)| = \J U > — u 2 -, an d 

N a 2 (a* O - a*) 
<3m X zz {Q^) = - V- 9{ U> -u) 6{u- U< ) . (38) 



u\ — u 2 



This susceptibility is very similar to that proposed in Ref. JT3| phenomeno 



logically, and also found for the long ranged spin-^ chain |14], pl| in Ref. [|T| 



with one important difference. The spectral weight here is dominatd by the 
upper threshold of the two parameter continuum m>, whereas the weight is 



peaked at the lower threshold u < in Ref. fL3 |. It is straightforward to see 
that if we choose 

2 ( n \ |sin(Q/2-g)| 

ot (q,Q-q) = v — , (39) 

Vsmg ysin(<5 - g) 
then on using Eq. (|37D, the weight is shifted to the bottom, and we get 

Qmx 22 (Q,u) = — 1 = 9( U> -u) e{u-u<) . (40) 
c - lu 2 — u 2 



"< 

With this choice, the static correlation function can be evaluated from Eq. 
(H|). We find 

leading to the asymptotic behaviour ~ (—1)™/ n at long distances. Indeed one 
can use the two parameters c and v in Eqs. PO} - pi) together with the various 
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sum rules known, in order to obtain very realistic structure functions which 
mimic the behaviour of the nearest neighbour Heisenberg antiferromagnet. 
At finite temperatures, we find from Eq. (|3~l"D in the usual case of a = 7 = 1 

< Sffi >= \6 r , - i [f n {^)f , (42) 

with 

fn(—) — — [ dx sin(nx) tanh(— sinx) , (43) 
2 n Jo 2 

leading to an exponentially decaying correlation function with a correlation 
length £ ~ 1/T for T — > 0. The function f n vanishes for even n in con- 
trast to one's usual expectation. In the presence of the phenomenological 
a, one must necessarily cut off the linear divergence of a at Q = it and 
q ~ 0, it. A temperature dependent cutoff, such as a 2 (a, b) = (| sin(a — b)/2\ + 
(const) 2 T)/( ^Jsm(a) + (const)Ty / 'sin(6) + (const)T ) interpolates nicely be- 
tween the zero temperature limit and the high temperature limit, and again 
gives a correlation length ~ 1/T. 

VI. Magnetic Fields 

We will now discuss the H-F ground state of the spin chain in the presence 
of uniform and staggered magnetic fields, and calculate the two susceptibili- 
ties. 

A. Uniform Magnetic Field 

For an uniform magnetic field Bz, we add a term —gi^BBJ2 n S'n to the 
Hamiltonian (|5|). Since this term commutes with (|5|), we can use the same 
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H-F decomposition as in fllQ| ) with g = 2 /it. Since the extra term in the 
Hamiltonian is quadratic in the Majorana operators, we only have to perform 
a rediagonalization of fllTl). We find that modes with S z = ±1 have an energy 

4J 

v±(q) = — sin g =f gi/i B B , (44) 

TT 

while the energy of the S z = modes remain unchanged. For B > 0, let us 
define a momentum q such that 

. -i ,7rgifj, B B 
1o = sin ( 4J ) , (45) 

and < q a < tt/2. (Such a q Q exists only if the magnetic field is less than a 
critical value B c = AJ/irgiyLB). Then the modes with S z = 1 and momenta 
lying in the range < q < q and tt — q Q < q < tt have negative energy, and 
the ground state of the system is one in which those modes are occupied. 
The change in the ground state energy is therefore given by a sum over all 
the occupied modes q, 

/ 4J \ 
AE Q HF = { — sin g - g x \i B B J 

q n 

AN J N 9U i B B 
= — 5- (1 - cosg D ) - q Q . (46) 

7T Z IT 

The expectation value of S z in the ground state is obtained either by counting 
the number of occupied modes, or by differentiating fl46|) with respect to 
gifi B B. Thus 

(S z ) = ^ = * sin-f^^). (47) 

7T TT \ 4J ' 

Finally, the (uniform) susceptibility is given by 

I = = ». (48) 



This agrees with the result in the previous section. For a strong magnetic 
field B > B c , the ground state is fully polarized with S z = N/2. These results 
are to be compared with the exact results for the susceptibility x — N/n 2 J, 
and the critical field B c = 1J j gi\x B |Tj|. 

Since S* has a non-zero expectation value in the ground state, the above 
calculation is not entirely self-consistent, i.e., one should also allow H-F de- 
compositions of the form 

(ft ft) = ifo 
and (M±i> = i/±i- (49) 

Further, the expectation values 

( K C +1 } = ( 4>i (Pl+i > = igr , 

and ( <p z n (f) z n+1 ) = ig L (50) 

may be unequal since the magnetic field breaks rotational invariance. On 
doing this more general H-F calculation, we find that although the ground 
state remains the same qualitatively (i.e., a number of S z = 1 modes have 
to be filled in the regions < q < q Q and tt — q Q < q < n), various numbers 
change. For instance, q Q is now given by 

q a + sin q a (1 + cos q ) = — — — . (51) 

The H-F parameters are 

2 2 

g T = cos q Q , g L = - , 

7T 7T 

f - 



f±i = . (52) 

7T 
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Since the magnetization is equal to Nq /n, the susceptibilty is x — N/6 J. 
(The critical field for complete polarization is B c = J(l + 2/n) / gin B ). We 
therefore have the curious result that a completely self-consistent H-F calcu- 
lation does not agree with linear response theory for small fields. 

B. Staggered Magnetic Field 

We now study the situation with a staggered magnetic field. We add 
a term —giiiBB^ n {—l) n S^ to the Hamiltonian and perform a H-F decom- 
position. As in the uniform case, we will assume that gx = gi = 2/7T 
and f = f±i = in Eqs. ([|9}|5C]) even though this is not completely self- 
consistent. We then find that the dispersion of the longitudinal modes remain 
the same as before while those of the transverse modes change. To be explicit, 

= — sing, 

7T 

16t7^ 1/2 
and u T (q) = ( —5- sin 2 q + gfn B B 2 ) . (53) 

Further, the change in the ground state energy is 

/ 47 \ 
AE HF = { — sin g - u T (q) J . (54) 

0<g<7r 71 

On differentiating this with respect to gifisB, we find the staggered magne- 
tization to be 

( E ("I)" Sn ) = NMbB /* g -1 (55) 
V Jo 2ir u T {q) 

For small fields, this goes as (NgifiBB /4J)ln(J / gifisB) which implies that 

the staggered susceptibility is divergent. This is the correct result. For large 

fields, the staggered magnetization approaches N/2 as it should. 
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VII. Discussion 



To summarize, we have used a Majorana fermion representation to study 
a nearest-neighbor isotropic antiferromagnetic spin-| chain. Within a trans- 
lation invariant Hartree-Fock approximation, we have found the spectrum of 
low-lying excitations, the two-spin correlation function, the structure func- 
tion, and the magnetic susceptibilities. All of these agree qualitatively with 
the results found earlier by a variety of other methods. The agreement can be 
made quantitative if we introduce some phenomenological functions within 
the Majorana formalism. 

It is somewhat surprising that a fully dimerized Hartree-Fock approxima- 
tion leads to a ground state with a lower energy. One way of stabilizing the 
translation invariant ground state with respect to the dimerized one is to ap- 
ply an uniform magnetic field with a strength B > 0.5829-B c = 0.7422 J/gipLs- 
Such a magnetic field lowers the energy of the translation invariant ground 
state below —3NJ/8, and does not change the energy of the dimerized ground 
state, for B < J/giUs, due to the finite gap to spin excitations. 

It would be interesting to go beyond our Hartree-Fock treatment and 
study the effects of fluctuations. Besides producing more accurate numbers 
for various quantities such as the spin wave velocity, such a study could also 
lead to a more detailed understanding of the "spinons" in a spin-| chain in 
terms of Majorana fermions. 

It may be instructive to examine models with anisotropy, frustration, and 
higher dimensionality using the Majorana representation, and to compare 
with known results. Amongst other things, this would help to determine the 
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range of validity of this way of studying spin-- systems. 

We have briefly examined the ferromagnetic case in which the exchange 
constant in Eq. (]5|) is negative. We perform a non-rotation invariant Hartree- 
Fock decomposition by allowing = —%<^ x n ^) v n to take an expectation value. 
We then obtain the correct ground state energy E Q = NJ/4, with the total 
S z = ±N/2. However we get the wrong dispersion relation, including a gap, 
for the low-energy excitations. Thus the Majorana Hartree-Fock approxima- 
tion is not a good starting point for studying the spin-^ ferromagnet. 
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